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LONG  TERM  GOALS 

The  development  and  improvement  of  cloud  microphysical  and  radiative  parameterizations  for  use  in 
mesoscale  models 

OBJECTIVES 

Investigation  of  marine  stratocumulus  clouds  microphysics  and  radiative  processes  using  the  CIMMS 
LES  model  with  explicit  microphysics  and  radiation.  The  data  from  FIRE  II/ASTEX  and  MAST  field 
experiments  will  be  used  to  validate  the  model  and  to  improve  our  understanding  of  the  interactions 
between  the  microphysical,  radiative,  and  thermodynamical  processes. 

APPROACH 

The  modeling  part  of  the  research  is  based  on  the  CIMMS  3-D  LES  model  of  boundary  layer 
stratocumulus  clouds  with  explicit  formulation  of  aerosol  and  cloud  drop  size-resolving 
microphysics.  The  model  has  been  thoroughly  tested  against  observations  from  ASTEX  and  MAST 
field  programs.  Use  the  model  to  generate  3D  data  fields,  including  the  rates  of  various 
microphysical  processes  needed  to  find  relations  between  bulk  variables  that  can  be  forecasted  in 
numerical  weather  prediction  models.  Use  the  3D  drop  size  distributions  to  calculate  the  optical 
properties  of  clouds  and  study  radiative  transfer  in  inhomogeneous  cloud  media  based  on  a  3D 
Monte  Carlo  radiative  transfer  model. 

WORK  COMPLETED 

A  parameterization  of  the  cloud  drop  activation  process  in  marine  stratocumulus  in  terms  of 
mesoscale  model  prognostic  variables  has  been  developed.  The  parameterization  accounts  for  the 
dependence  on  CCN  concentration,  as  well  as  intensity  of  turbulence  in  the  boundary  layer. 

A  design  of  a  new  parameterization  of  cloud  physics  processes  based  on  the  integral  moments  of  the 
drop  spectra,  as  opposite  to  the  partial  moments  used  in  Kessler-type  parameterizations  have  been 
developed.  The  importance  of  probability  distribution  functions  (pdf)  of  prognostic  variables  in 
applying  the  cloud  physics  parameterization  in  heterogeneous  broken  cloud  fields  has  been 
demonstrated  using  LES  model  data. 

The  effect  of  inhomogeneous  cloud  media  on  3D  radiative  transfer  was  studied  using  the  coupled  3D 
Monte-Carlo  radiative  transfer  and  explicit  microphysics  cloud  models.  We  have  showed  the  impact 
of  cloud  vertical  inhomogeneity  on  satellite  retrieval  of  cloud  fraction. 
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RESULTS 


a.  Parameterization  of  the  cloud  drop  activation  process 

The  number  of  activated  cloud  drops,  NCd  is  parameterized  in  terms  of  the  domain  average  number  of 
CCNs  at  1%  supersaturation,  Nccn,  and  turbulence  intensity  characterized  by  the  domain  average 
vertical  velocity  variance,  <w2>.  We  have  considered  environments  Nccn  in  the  range  from  25  to  1000 
cm"3  and  values  of  <w2>  in  the  range  from  0.06  to  0.27  m2/s2.  The  drop  concentration  averaged  over 
the  whole  cloud  layer  was  related  to  Nccn  and  <w2>  using  regression  analysis,  resulting  in  the 
following  link: 


Ncd  =N™2  <w2> °-15 

Fig.  1  shows  the  number  of  activated  drops  for  three  levels  of  turbulence  intensity  depicted  by 
different  symbols.  The  best  fit  to  the  empirical  data  obtained  in  various  field  programs  by  Martin  et  al 
(1994)  and  O’Dowd  et  al  (1996)  is  also  shown.  The  derived  parameterization  agrees  well  with  both 
empirical  fits  in  the  range  of  CCN  concentrations  up  to  500  cm"3.  The  latter  covers  the  most 
commonly  observed  marine  environments.  The  dependence  of  drop  concentration  on  the  intensity  of 
vertical  motions,  <w2>,  is  stronger  for  more  polluted  environments.  The  mesoscale  grid  averaged 
ambient  aerosol  concentration,  as  well  as  the  level  of  turbulence  intensity,  is  available  in  the  current 
version  of  the  ARPS  model  where  the  parameterization  has  been  already  implemented.  The 
parameterization  can  also  be  easily  included  in  the  Navy  COAMPS  model. 
h.  Effect  of  sub  grid  scale  inhomogeneity 

The  conversion  rates  that  form  the  basis  of  a  cloud  physics  parameterization  are  derived  and  applied 
locally,  i.e.,  relating  variables  defined  on  a  cloud  scale.  The  extrapolation  of  a  cloud-scale 
parameterization  to  a  mesoscale  model  requires  empirical  knowledge  or  theoretical  prediction  of 
statistical  distributions  of  cloud  prognostic  variables,  e.g.,  pdf.  The  scale  dependence  of  cloud 
parameterization,  as  well  as  the  effect  of  partial  cloud  cover  in  a  grid  box  has  been  emphasized  in 
many  studies  (e.g.,  Sundqvist  1993;  Randall  1995).  Even  for  cases  with  total  cloud  cover,  the  neglect 
of  cloud  micro  structure  inhomogeneity  may  lead  to  serious  errors  in  grid  average  values. 

We  demonstrated  the  impact  of  cloud  microstructure  inhomogeneity  by  expanding  the  cloud-scale 
stratocumulus  parameterization  of  Khairoutdinov  and  Kogan  (1998a)  to  a  mesoscale  grid  of  several 
kilometers.  As  an  example  we  used  the  autoconversion  rate  from  the  above-mentioned 
parameterization  which  is  a  non-linear  function  of  cloud  prognostic  variables. 

The  autoconversion  rate  averaged  over  the  mesoscale  grid  is  given  by 

^fYaut--, /  9(2  e,Nc)  f(Qc,Nc)  dQJN, 

where  cp(g(;,AQ  is  the  joint  pdf  of  Q,  and  N,  and  f(Qc.Nc )  is  the  cloud-scale  autoconversion  rate. 

Table  1  shows  the  average  autoconversion  rate  calculated  using:  1)  the  joint  pdf(p(Qc,Nc ) ;  2)  the 
separate  pdf  for  Q,  and  Nc\  and  3)  the  domain  averaged  values  <QC>  and  <NC>.  As  one  might  expect, 
the  autoconversion  rate  depends  significantly  on  the  type  of  pdf  used  in  calculations  (joint  pdf  versus 
separate  pdf)  and  varies  with  time  reflecting  the  evolution  of  drizzle  in  the  simulation  of  a 
stratocumulus  cloud  layer  and  the  corresponding  changes  in  pdf . 


Table  1.  The  dependence  of  autoconversion  rates  on  subgrid  statistics 


Time,  hrs 

Joint  pdf 

Separate  pdf  for  Q  &  N 

Constant  <Q>  and  <N> 

3.0 

2.88x10  s 

3.72xl0"8 

0.61x10  s 

4.5 

8.82xl0"s 

ll.llxl  0”8 

0. 58x10  s 

6.0 

21.25xl0"8 

13.88x10  s 

0.90x10  s 

The  above  example  emphasizes  the  need  for  experimental  studies  of  statistical  distribution  of  cloud 
parameters  inside  a  mesoscale  grid  box.  Another  possibility  to  obtain  the  subgrid  statistics  is  to  use 
the  grid  embedded  models  (GEM).  The  rapid  progress  in  computer  power  makes  practical  the 
incorporation  of,  at  least,  a  2D  cloud-resolving  model  into  a  fraction  of  a  mesoscale  model  grid  cells. 
The  resolvable  variables  of  the  mesoscale  model  can  serve  as  the  ambient  parameters  for  the  GEM, 
while  the  pdf  and  other  subgrid  characteristics  calculated  by  the  GEM  (e.g.,  surface  fluxes,  TKE)  will 
enter  the  parameterizations  used  in  the  parent  mesoscale  model. 

c.  Investigation  of  the  role  of  sub  grid  inhomogeneity  on  radiation  parameterizations 
This  study  explores  how  the  grid  averaged  radiative  parameters  depend  on  cloud  micro  structure 
vertical  inhomogeneity.  The  radiative  parameters  are  calculated  by  the  3D  Monte  Carlo  radiative 
transfer  model  that  uses  as  input  the  3D  field  of  cloud  drop  distributions  provided  by  the  LES  model. 
The  difference  between  the  fully  inhomogeneous  benchmark  case,  3D,  and  the  vertically 
homogeneous  case,  2D,  (see  Fig.  2)  illustrates  the  problem  of  cloud  satellite  retrieval  of  cloud 
fraction.  The  retrieval  algorithms  classify  pixels  as  cloudy  or  non-cloudy  depending  on  the  threshold 
value  of  albedo  Rc.  The  pixels  with  albedo  exceeding  Rc  are  flagged  as  cloudy,  while  those  with 
albedo  less  than  Rc  as  clear-sky.  As  evident  from  Fig.  2,  there  is  a  much  larger  number  of  non-cloudy 
pixels  in  the  C2D  case  than  in  the  C3D  case.  This  suggests  that  neglect  of  vertical  stratification  of  the 
extinction  coefficient  may  lead  to  overestimation  of  cloud  fraction.  Our  results  also  show  that 
vertical  profiles  of  both  the  mean  and  the  variance  of  the  extinction  coefficient  are  important  in  cloud 
fraction  retrieval. 

IMPACT 

The  improved  parameterization  of  the  physical  processes  in  marine  stratocumulus  clouds  will  result 
in  more  accurate  weather  prediction  for  Navy  operations.  In  particular,  the  work  is  aimed  at 
improved  prediction  of  atmospheric  visibility,  precipitating  cloud  layers,  and  cloud  optical  and 
radiative  parameters. 

TRANSITIONS 

Our  work  is  known  to  the  COAMPS  development  team  at  NRL.  The  developed  parameterizations 
are  planned  to  be  implemented  in  the  COAMPS  model  as  part  of  the  MURI  Grant  to  the  University 
of  Oklahoma. 

RELATED  PROJECTS 

The  current  proposal  is  aimed  at  development  of  physical  parameterizations  for  cloud  scale  (LES) 
models.  It  is  closely  related  to  the  ONR  project  “Remote  sensing  and  prediction  of  the  coastal 
marine  boundary  layer”  (N00014-96-1-1112).  The  latter  project’s  goal  is  to  formulate  the 
parameterizations  for  application  to  mesoscale  prediction  models. 


REFERENCES 


Khairoutdinov,  M.  and  Y.  L.  Kogan,  1998a:  A  new  cloud  physics  parameterization  for  large-eddy 
simulation  models  of  marine  stratocumulus.  Mon.  Wea.  Rev.,  conditionally  accepted. 

Martin,  G.  M.,  D.  W.  Johnson,  and  A.  Spice,  1994:  The  measurment  and  parameterization  of 
effective  radius  of  droplets  in  warm  stratocumulus  clouds.  J.  Atmos.  Sci.,  51,1823-1842. 

O'Dowd,  C.  D.,  M.  H.  Smith,  I.  E.  Consterdine,  and  J.  Lowe,  1997:  Marine  aerosol,  sea-salt,  and  the 
marine  sulphur  cycle:  a  short  review.  Atmos.  Env.,  73-80. 

Randall,  D.  A.,  Parameterizing  fractional  cloudness  produced  by  cumulus  detrainment.  Paper 
presented  at  Workshop  on  cloud  microphysics  parameterizations  in  global  atmospheric  circulation 
models.  23-25  May  1995,  Kananaskis,  Alberta,  Canada.  WMO/TD-No.  713,  1-16. 

Sundqvist,  H.,  1993:  Parameterization  of  cloud  in  large-scale  numerical  models.  In  Aerosol-Cloud- 
Climate  Interactions.  P.  V.  Hobbs  (ed.),  Academic  Press.  175-203. 

PUBLICATIONS 

Ferek,  R.  J.,  T.  Garrett,  S.  Strader,  K.  Nielsen,  G.  E.  Innis,  J.  P.  Taylor,  A.  S.  Ackerman,  Y. 

Kogan,  Q.  Liu,  B.  A.  Albrecht,  D.  Babb,  1998:  Drizzle  Suppression  in  Ship  Tracks.  J.  Atmos.  Sci., 
accepted. 

Liu,  Q.,  Y.  L.  Kogan,  D.  K.  Lilly,  M.  P.  Khairoutdinov,  1997:  Variational  optimization  method  for 
calculation  of  cloud  drop  growth  in  an  Eulerian  drop-size  framework,  J.  Atmos.  Sci.,  54,  2493-2504. 
Khairoutdinov,  M.  F.  and  Y.  L.  Kogan,  1998:  A  Large  Eddy  Simulation  Model  with  Explicit 
Microphysics:  Validation  Against  Aircraft  Observations  of  a  Stratocumulus-Topped  Boundary 
Layer.  J.  Atmos.  Sci.  accepted. 

Khairoutdinov,  M.  and  Y.  L.  Kogan,  1998:  A  new  cloud  physics  parameterization  for  large-eddy 
simulation  models  of  marine  stratocumulus.  Mon.  Wea.  Rev,  conditionally  accepted. 

Liu,  Q.,  Y.L.  Kogan,  D.K.  Lilly,  D.W.  Johnson,  G.E.  Innis,  P.A.  Durkee,  K.  Nielson,  1998:  LES 
modeling  of  ship  track  formation  and  its  sensitivity  to  boundary  layer  structure.  J.  Atmos.  Sci. , 
accepted. 

IN-HOUSE/OUT-OF  HOUSE  RATIOS 


100%  of  work  is  done  by  out-of-house  organization  (University  of  Oklahoma) 


00 

I 

S 

CJ 


C 

o 

cd 

-+j 

CD 

a 

£ 

o 

o 

Qh 

o 

S-H 

Td 

td 

2 

o 

G 


Fig.  1.  The  concentration  of  activated  cloud  drops  as  a  function  of  CCN  concentration  and 
turbulence  intensity.  The  LES  parameterization  is  shown  by  the  dashed  line,  three  levels  of 
turbulence  intensity  are  depicted  by  different  symbols.  The  best  fit  to  the  empirical  data  obtained 
from  Martin  et  al  (1994)  and  O’Dowd  et  al  (1996)  is  also  shown. 
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Fig.  2.  Spatial  distribution  of  the  albedo  field  for  the  inhomogeneous  cloud  (right,  3D)  and  the 

vertically  homogeneous  cloud  (left,  2D). 


